Hidden Scaling in the Quantum Hall Metal-Insulator Transition 
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Scaling properties of the quantum Hall metal-insulator transition are severely affected by finite size 
effects in small systems. Surprisingly, despite the narrow spatial range where probability structure 
functions exhibit multifractal scaling, we clearly verify the existence of extended self-similarity - a 
hidden infrared scaling phenomenon related to the peculiar form of the crossover at the onset of 
nonmultifractal behavior. As finite size effects get stronger for structure functions with negative 
orders, the parabolic approximation for the multifractal spectrum loses accuracy. However, by 
means of an extended self-similarity analysis, an improved evaluation of the multifractal exponents 
is attained for negative orders too, rendering them consistent with previous results, which rely on 
computations performed for considerably larger systems. 
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The quantum Hall metal-insulator transition has been 
the focus of great interest, proving to be an ideal stage for 
the interplay of several approaches in the physics of con- 
densed matter localizationi^ In fact, while a number of 
field theory descriptions have been proposed ,^'*^'^ numer- 
ical simulations and experiments are relatively precise;^ 
and reduced lattice modelsSii amenable to theoretical 
treatment, have an interesting connection with the clas- 
sification problem of disordered systems.' Even though 
the basic physical picture of the quantum Hall metal- 
insulator transition was put forward around two decades 
ago,^ a lot of work has been devoted to the subject 
along the last years. In particular, the complete the- 
oretical characterization of the multifractal behavior at 
the transition,^*' and an exact solution, likely in terms of 
some conformal field theory,- are so far open problems. 

Our aim in this work is to call attention to a hid- 
den scaling property of the non-interacting quantum 
Hall metal-insulator transition, named extended self- 
similarity (ESS). This kind of scaling phenomenon was 
found for the first time in the context of fully devel- 
oped turbulence^^ii^ where it has been a fundamen- 
tal tool to extract information on the turbulent regime 
out of moderate Reynolds number data. After the 
pioneering works of Benzi et al., ESS has been veri- 
fied in other classical multifractal dynamical systems, 
as chaotic advection^, kinetic roughening, turbulent 
magnetohydrodynamics^ and strange attractors.^^ Re- 
garding the quantum domain, a clear numerical evidence 
of ESS in the two-dimensional metal-insulator transi- 
tion has been recently reported, -'^^ for a model of random 
Dirac fermions, which is a close variant of the quantum 
Hall problem^i^ 

In short words, the ESS phenomenon corresponds to a 
specific form of crossover for the scaling functions, that 
occur between small and large length scales. More con- 
cretely, let L be the linear system size, ipi^) be the wave- 
function for some fixed energy level and O be a square 
with size r x r. We define the probability structure func- 



tion of order q as 



(1) 



where the average is taken over disorder realizations. The 
multifractal scaling law Sq{r) ~ r''^'^^ supposedly holds at 
the metal-insulator transition. However, if finite size ef- 
fects are taken into account or if the disordered potential 
has a finite correlation length, the scaling law is valid 
only through a certain range of scales. Either towards 
the smaller or larger length scales (the "ultraviolet" and 
"infrared regimes", respectively), wavefunctions will be- 
come effectively smoothjil leading to the trivial scaling 
({q) = 2q. As r moves out from the scaling region, to- 
wards the infrared, the general meaning of ESS is just a 
peculiar profile of scaling, in the formic 



5,(r)=c,[/i(r)]C(^)[/2(r)] 



2q 



(2) 



crossover functions introduced above behave as 
= r/L and /2(?') = 1 in the scaling region and 
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as /i(r) = 1 and f2ir) ~ r/L in the infrared regime, 
where finite size effects become relevant. 

In the particular problem of localization, it follows 
from translation invariance and from definition 1^ that 
5i(r) = {r/L)^ is an exact relation. Taking this result 
into account together with Eq. we find that 



Cg[/l(0]'^'^^"'^ 



(3) 



Therefore, even if a plot of ln[5g(r)] against ln(r) does 
not exhibit scaling when larger length scales are included, 
a plot of \ia[Sp{r)r~^P] against \n[Sq{r)r~^'^] will do: nu- 
merical data collapses on a straight line, with slope 



dp) - 2p 
Ciq) - 2q 



(4) 



As a consequence, experimental or numerical evaluations 
of ri{p, q) would be subject to smaller errors than the 
scaling exponents of probability structure functions. 



2 



In order to investigate the issue of ESS, we model the 
quantum Hall metal-insulator transition in a square lat- 
tice through the nearest-neighbor hamiltonian 

H = J2 ^^4c^ + E ^(e'*'^4c, + e-'*--4c.) . (5) 

Defining the lattice parameter as the unit length, the 
Landau gauge gives (jjij — along the x direction and 
= ±27ra;'i>/$o along the =Fj/ direction, where the mag- 
netic flux ratio is $/$o = Be/h. On-site disorder is 
provided by the random energies Si, correlated within a 
length scale of order A. In practice, they may be de- 
fined as Si — I/ttA^ exp(— jxi — S'jl^/A^), in terms 
of uncorrelated random variables a^, with |ai| < 11^/2. 

The simulations were taken for a 60 x 60 lattice, with 
periodic boundary conditions and $/$o = 1/20. For this 
choice of parameters, the linear dimension of the system 
is L ~ 34€, where £ — y/h/eB is the magnetic length. 
Furthermore, we set A = 2, the hopping parameter V — 
0.14, and W/V = 1.6. 
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FIG. 1: The averaged participation ratio as a function of 
energy for the first Landau level. The averages are taken over 
45 disorder realizations in a 60 x 60 lattice. The system's 
linear size is L ~ 34£. 

According to the generally accepted picture, the large 
quantum degeneracy of Landau levels is broken due to 
disorder, and delocalized states are found only around a 
critical energy close to each band center. A careful deter- 
mination of the critical energy, which takes into account 
the levitation of extended states f2£ is based on the com- 
putation of the averaged participation ratic^^'^^ (PR); 
written in terms of the wavefunction amplitude as 

In our computations of the participation ratio, which 
are restricted to the first Landau level and shown in Fig.l, 
the total number of sites is = 60 x 60, and the above 
average is taken over a set of 45 disorder realizations. 
For each one of them, the hamiltonian Q was exactly 
diagonalized, producing 180 (= 60 x 60 x 1/20) states 



per Landau level. The more localized is a state the lesser 
is its participation ratio. 

Thus, to study the statistical properties of the delocal- 
ized wave functions, it is necessary to perform a selection 
of states around the participation ratio peak. We focused 
on the energy range —3.703 < E/V < —3.695, which 
contains 8 states. We picked up, out of these states, the 
one with the maximum PR for each disorder realization. 
Had we chosen all of them for each realization, the sta- 
tistical analysis would be spoiled by the several localized 
states that eventually appear, as a consequence of the 
strong fluctuations inherent to small systems. In fact, 
their presence can be inferred from a comparison of dif- 
ferent selections. The selected PR's fall in the interval 
0.34 ± 0.04. An alternative (but useless) selection which 
takes the minimum PR for each realization gives PR's 
in the interval 0.09 ± 0.06. On the other hand, taking 
all of the 8 states we get 0.23 ± 0.1. Of course, these 
three types of selection will lead to converging results in 
the large box limit L ^ oo. We assume here that our 
prescription is the one which yields the fastest conver- 
gence - an hypothesis we verified a posteriori, through 
the computation of the multifractal exponents Cil)- 
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FIG. 2: The structure function of order q=10 reveals, when 
multiplied by r"^" (circles), a scaling region (6 < r < 14) 
fit by a straight line with slope —8.77. The crossover to a 
less fluctuating regime (14 < r < 31) is observed, with slope 
— 5.72. The inset shows the ultraviolet, scaling and the in- 
frared regions, covering the interval 2 < r < 51. The fltting 
in the inset corresponds to the steeper straight line in the 
large plot. 

It is interesting to observe that the energy position 
of the participation ratio peak reaches in a fast way its 
asymptotic valued showing little additional levitation 
when the system's size becomes larger than ^ 40 x 40, 
while the width and height of the PR curve are still sen- 
sitive to finite size effects. 

We were able to find the scaling region of probability 
structure functions through an analysis of S q{r)r~''^'^ , as 
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shown in Fig. 2. In the uhraviolet and in the infrared 
scales, as aheady discussed, this function becomes flatter. 
The scaling region, therefore, appears as an inflection 
between these two regimes. This fact is shown in the 
inset of Fig. 2 for the structure function of order q — 10. 
In general, for positive orders, scaling takes place for 6 < 
r < 14 (r is measured in lattice units). A crossover to the 
infrared regime is observed for 14 < r < 31. In Fig. 2, we 
establish fittings with slopes -8.77 (scaling region) and - 
5.72 (larger length scales). These numbers have a relative 
variation of ~ 53%. Similar large variations were found 
for other positive orders, in the interval < g < 14. 
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FIG. 3: Structure functions of different orders are compared 
and the phenomenon of extended self-similarity is clearly ob- 
served. The fitting, which has slope r;(10,3) = 8.37, approxi- 
mates the whole numerical curve, and is produced taking only 
the data for 6 < r < 14 (i.e., the first 9 circles from right to 
left). 

The existence of ESS is strikingly indicated in Fig. 3, 
the central result of our work, where we compare struc- 
ture functions with orders p — 10 and q — 3, to obtain 
ry(10, 3). The straight line corresponds to the fitting for 
6 < r < 14, which has slope 8.37. The fitting at larger 
length scales, for 14 < r < 31, has slope 8.50. Here, the 
slopes have a relative 1.5% variation, which should be 
compared to the 53% variation obtained in Fig. 2. Such 
higher precision evaluation is the hallmark of ESS, which 
was also verified for other pairs of orders (p, q) . A rel- 
evant question is how small or how disordered can be 
the system so that ESS is still meaningful. While we 
do not have yet a conclusive answer to this point, we 
have inspected simulations on a 40 x 40 lattice, getting 
results which are essentially equivalent to the ones de- 
scribed here. 

We have checked that ESS keeps holding for neg- 
ative orders, where however the parabolic approxima- 
tion for the multifractal spectrum becomes inaccurate 
(as opposed to positive orders). The parabolic ap- 



proximation simply consists of assigning a log-normal 
probability distribution function Pr(ck) ~ r^~f(°')^ with 
f{a) = 2 — (a — ao)^/4(ao — 2), to have a box probability 

(Px\tlj{x)\'^ that scales as r", for a given a parameter 
The function f{a) is the so-called multifractal spectrum, 
and can be interpreted as the fractal dimension of an "iso- 
a set" . In this framework, the multifractal exponents of 
probability structure functions will depend solely on a 
single parameter ao. They are given by 

aq)=2q + qiao-2)il-q) . (6) 

This quadratic expression is assumed to hold as long 
as the multifractal spectrum is positive at the saddle- 
point of r'^°'pr{a), i.e., we need to have /(a) > 0, 
with a = ao — 2q(aQ — 2). Numerical simulations, in- 
volving different types of disorder in systems with lin- 
ear sizes L ~ 200£, yield ao = 2.25 ± 0.05 for the 
first Landau level. More recent numerical work, per- 
formed in the Chalker-Coddington model for lattices with 
sizes up to 1280 x 1280, have advanced the conjecture 
that the parabolic approximation is in fact exact, with 
ao = 2.261 ± 0.003,23: 

As it can be easily checked, it is worth noting that in 
the parabolic approximation, we have 



which is remarkably independent on uq. Computations of 
7y(p, q) could be used as a rigorous test of the parabolic 
approximation. However, it is important to point out 
that small errors in the multifractal exponents may lead 
to relatively larger errors for rj^p, q). 

We show in Fig. 4 the multifractal exponents ({q) for 
—2.5 < q < 2.5, used to analyze the parabolic approxi- 
mation for the multifractal spectrum. Circles represent 
the direct numerical evaluations. The spatial range where 
ESS is now verified, viz 6 < r < 16, is smaller compared 
to the one related to positive orders. The exponents plot- 
ted in Fig. 4 (circles) correspond to linear fittings in the 
interval 6 < r < 10. If the positive orders are neglected, 
the best fit produced from Eq. © would give ao = 2.85, 
which lacks satisfactory precision. On the other hand, if 
the negative orders are neglected, the best fit is found for 
ao = 2.2, now a tolerable value. We conclude in this way 
that stronger finite size effects occur for the case q < 0. 

There is in fact a simple phenomcnological interpreta- 
tion of the above result, along the lines of the cascade 
picture of localization addressed in Ref. [17]. There, 
a random /3-model approach has been implemented, al- 
lowing us to write the multifractal scaling exponents as 
C{q) — 2q ~ \og^{(3~'^). The random parameter < (3 < 
oo is used to model regions of expanding/collapsing evo- 
lution of the wavefunction support, in a discrete iteration 
defined by the scale compression parameter a > 1 (i.e., 
sub-regions with linear sizes Lo, Lo/a, Lq/o^, ... are 
successively generated). Essentially, the iteration of a 
given "coherent" region which belongs to the wavefunc- 
tion support is identified to a mapping with local jacobian 
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/3. Therefore, taking g < in the latter expression for 
C{q), it is clear that (3 > 1 becomes relevant - and so do 
expanding fluctuations of the wavefunction - which we 
expect will lead to stronger finite size effects. 
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FIG. 4: The numerical scaling exponents for the probabil- 
ity structure functions (circles), with orders —2.5 < q < 2.5. 
The solid line corresponds to the theoretical exponents pre- 
dicted from the parabolic approximation for the multifractal 
spectrum, with ao = 2.2. The squares are the estimated cor- 
rect values of exponents for negative orders, taking ESS into 
account. 

Even though the values of C('z) for g < strongly de- 
viate from the parabolic approximation predictions, we 
have verified that the ESS slopes are close to the exact 
values: taking, for instance, p = —2.5 and q = —1.5, 
we get rj{—2.5, —1.5) = 2.25, just 4% off the exact value, 
obtained from (TJ, which is 7/3 ~ 2.33. This suggests 
that the numerical results for ({q) — 2q are proportional 



to the exact ones. Taking this hypothesis into account, 
we have devised an optimization procedure whereby the 
solid line (obtained from the parabolic approximation) 
and the squares (a recovery of more precise multifractal 
exponents) are found simultaneously in Fig. 4. We replace 
the multifractal exponents for q < from their measured 
values Ciq) to c(C(g) — 2q) + 2q. The parameters c and 
are chosen to be the ones that yield the minimum devia- 
tion from the solid line. We find in this way c = 0.23 and 
ao = 2.2. If a more restrict range of g's is considered, 
say —1.5 < q < 1.5, we find c = 0.36 and ao = 2.27. 
Thus, we estimate ao = 2.24 ± 0.04 which is in excellent 
agreement with previous known results. Observe that the 
only bias in the method is the fact that the multifractal 
exponents for positive orders were preserved. 

In summary, we have numerically established the exis- 
tence of ESS in the quantum Hall metal-insulator transi- 
tion for a small system (L ~ 34£) , where strong finite size 
effects cannot be avoided. We observed that the scaling 
behavior of structure functions with negative orders is 
not derived from the parabolic approximation. However, 
an optimization scheme based on ESS considerations was 
introduced, which led to a more precise determination 
of the scaling exponents for negative orders. We found 
an agreement with known standard results, obtained on 
their turn from the analysis of much larger systems. We 
hypothesize that the ESS phenomenon could be used as 
a valuable tool in order to improve the scaling exponents 
presently reported in the literature for the quantum Hall 
effect. A theoretical foundation of ESS is highly desir- 
able, likely through field theory methods, which could 
allow to check if the hidden scaling behavior found in 
this work is indeed a general property of metal-insulator 
transitions. 

This work has been partially supported by CAPES, 
FAPERJ and FAPESP. 
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